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Thermalization of a parametrically driven magnon gas leading to the formation of a Bose-Einstein conden¬ 
sate at the bottom of a spin-wave spectrum was studied hy time- and wavevector-resolved Brillouin light scat¬ 
tering spectroscopy. It has been found that the condensation is preceded by the conversion of initially pumped 
magnons into a second group of frequency degenerated magnons, which appear due to parametrically stimulated 
scattering of the initial magnons to a short-wavelength spectral region. In contrast to the first magnon group, 
which wavevectors are orthogonal to the wavevectors of the magnons at the lowest energy states, the secondary 
magnons can effectively scatter to the bottom of the spectrum and condense there. 


Parallel parametric pumping Qi, where an alternating 
magnetic pumping field h(f) is applied parallel to the direc¬ 
tion of a bias magnetic field H, can effectively couple the spin 
system of a magnetic material with an external microwave en¬ 
ergy source, and thus is widely used to generate (3,01, am¬ 
plify lIsl-Bl, and restore ||8l0 spin-wave signals in macro- and 
micro-sized magnetic structures. In addition, this technique 
fruitfully serves for experimental studies of nonlinear dynam¬ 
ics in multi-mode wave systems as well. One of the great¬ 
est achievements in this area was the realization of a Bose- 
Einstein condensation of magnons i[ii, where the parametri¬ 
cally pumped magnon gas undergoes a phase transition lead¬ 
ing to the spontaneous formation of a coherent state at the bot¬ 
tom of the spin-wave spectrum iQ . Since the moment of this 
discovery, the physics of the parametrically driven magnon 
gases is under active theoretical and experimental investiga¬ 
tion. In particular, the time dependent behavior of the magnon 
gas in the phase-energy space 112 ■Q, the temperature of the 
magnongas lll6|] , and the role of different scattering mecha¬ 
nisms fllTI] in the process of the thermalization of the injected 
magnons are in the focus of attention. 

Here, we provide experimental insight into the evolution of 
a magnon gas affected by four-magnon scattering in the pres¬ 
ence of an external pumping field. We claim that the magnons 
initially pumped to the transversal branch of the spin-wave 
spectrum (magnon wavevector q _L H) do not scatter into 
the lowest energy states located at the longitudinal branch 
(q II H). Such a scattering and the consequent formation 
of a Bose-Einstein condensate (BEC) are associated with the 
parametric excitation of a second group of short-wavelength 
magnons propagating at an angle Z(q, H) < 90°. Simulta¬ 
neously, the initially excited transversal magnon group is sup¬ 
pressed to the near thermal level by a mutual action of the 
second magnon group and the external pumping field. 

The measurements were performed using a low-damping 
ferrimagnetic film of yttrium iron garnet (YIG, Y 3 Fe 50 i 2 ) by 
means of a combined microwave and Brillouin light scattering 
(BLS) setup schematically shown in Eig.fTJa). The YIG film 
of 6.7 fim thickness with saturation magnetization AttMs — 
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FIG. I. (color online), (a) Sketch of the experimental setup. H is 
a bias magnetic field, qy and are in-plane wavevectors of the 
transversal (q _L H) and longitudinal (q || H) magnons. The setup 
is tiltable by the angles 0^^ and 0q|| allowing for wavevector reso¬ 
lution. (b) The magnon spectrum calculated for the first 38 thickness 
modes n Q. The shadowed area between the red arrows illustrates 
the magnon injection by parametric pumping. The initially pumped 
magnon group is marked by two violet dots. The second magnon 
group is shown by four green dots. The green arrows indicate group 
velocities of the secondary magnons. The yellow arrows outline the 
parametrically stimulated magnon scattering. 


1750G was grown in the (111) crystallographic plane on a 
gallium gadolinium garnet substrate by liquid phase epitaxy. 
The in-plane magnetized YIG-film sample with dimensions 
2x10 mm^ was placed on top of a 50 fim wide microstrip res¬ 
onator, which was used to induce the pumping Oersted field 
h(f) [see Eig.llfa)]. The resonator was driven by microwave 
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FIG. 2. (color online). Dependence of the density of parametrically excited magnons at half of the pumping frequency /p/2 (red dots) and 
the density of the magnons at the bottom of the spin-wave spectrum (blue squares) on the pumping power Pp. The magnon densities are 
measured as the intensities of the corresponding BLS signals. The selected temporal profiles of these signals are presented in panels (a)-(f). 
The pumping pulse of 10 /is duration is applied at zero moment of time. The data marked by the red dots and by the filled blue squares were 
measured from the corresponding temporal profiles just before switching off the pump pulse. The open blue squares show the intensity of the 
peak of the magnon Bose-Einstein condensate, which is formed due to the evaporative supercooling Gl of the magnon gas after the pumping 
is shut down [see e.g. panels (e)-(f)]. The letters near the data points relate them to the BLS temporal profiles in the corresponding panels. 
Three waveforms in panel (a) demonstrate the development of the conventional parametric instability with increasing pump power. 


pulses at a carrier frequency of /p = ujp/2Tr = 13.74 GHz 
with peak powers Pp ranging from 1 mW to 4 W. The chosen 
pump pulse duration of 10 p,s was sufficiently long to observe 
the response of the spin-wave system even at the lowest ap¬ 
plied pumping powers. A pulse repetition time of 100 ^s ex¬ 
cluded microwave-heating effects. 


In course of the parallel pumping process the photons of 
the microwave pumping field split into magnon pairs with op¬ 
posite wavevectors at half of the pumping frequency as illus¬ 
trated in Fig.[Tfb). The strength of the bias field H = 1735 Oe 
has been chosen to allow for a pumping of the magnon pairs 
slightly above the ferromagnetic resonance frequency (FMR). 
In this case, firstly, the parallel pumping achieves its highest 
efficiency because the magnons are pumped to the transversal 
spectral branch, which has the lowest threshold of the para¬ 
metric instability S. Secondly, no kinetic instability pro¬ 
cess, which corresponds to a one-step scattering of the para¬ 
metric magnons to the bottom of the spin-wave spectrum 
and which may disturb the formation of the magnon BEC, is 
allowed due to energy and momentum conservation 112 111 . 


The magnon system was analyzed by means of time- 
resolved BLS spectroscopy 11221] with a time resolution of 1 ns. 
A probing laser beam with a power of 5 mW and a wave¬ 
length of 532 nm was focused to a focal spot with a diameter 
of 50 fj,m in the middle of the microstrip resonator. The op¬ 
tics to focus the BLS laser onto the YIG sample was also used 
to collect the scattered light, which was analyzed by a multi¬ 
pass tandem Fabry-Perot interferometer. The magnons at the 
frequency /p/2 and at the bottom of the spin-wave spectrum 
were detected simultaneously. The corresponding BLS sig¬ 
nals were collected in 150 MHz wide frequency bands. The 
wavenumber sensitivity in the range ±5-10^ rad/cm was in¬ 


sured by use of a wide-aperture optics. 

Figure|2]presents the dependencies of the magnon densities 
in the chosen spectral regions on the input pumping power. 
The red and blue waveforms show the time evolution of the 
parametric and the magnons at the bottom of the spectrum, 
respectively. One can see in Fig. |2la) that the parametric 
magnons appear already at a small pumping power of 4.2 mW. 
The increase of the pumping power to 6.7 mW and to 9.6 mW 
results in an earlier appearance and in a faster saturation of 
the parametric magnons at higher density levels as shown in 
Fig. |2|a). In accordance with the theory of the parametric 
spin-wave instability iEbIi this saturation is caused by a non¬ 
linear dephasing of the excited magnons relative to the mi¬ 
crowave pumping field h(f). 

Surprisingly, a further increase of the pumping power leads 
to a consecutive decrease of the saturated magnon density. 
At a pumping power Pp ~ 1270mW, when the thermal- 
ized magnons start to be detected at the bottom of the spin- 
wave spectrum, the parametric magnons are hardly visible 
anymore [see waveforms in Fig.|2l[d)]. For pumping powers 
above 2770 mW the parametric magnons are practically non- 
detectable with the exception of a sharp initial peak similar to 
one in Fig.|2|d). Nevertheless, at Pp > 2770 mW the density 
of the low-energy magnons sharply rises up due to the forma¬ 
tion of the Bose-Einstein condensate. 

Conventional four-magnon scattering processes, which 
were considered in the past to be responsible for the thermal- 
ization of a magnon gas and, thus, for the formation of the 
magnon BEC, cannot be assumed as a reason for the disap¬ 
pearance of the parametrically pumped magnon group. Dur¬ 
ing the action of the external pumping this nonlinear mecha¬ 
nism can only limit but not reduce the number of the paramet- 
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ric magnons: The decrease in their density would lead to the 
consequent decrease of the scattering efficiency and, thus, to 
the density saturation. 

The phenomenon can be understood if one takes into ac¬ 
count the role of the external pumping field in the four- 
magnon scattering process. The time-dependent evolution of 
a spin wave with amplitude Cq, wave vector q, and complex 
frequency ujq = ojq + iVg can be described by the follow¬ 
ing equation 12311 . where is the spin-wave relaxation fre¬ 
quency: 


dc, 

dt 


=iU}qCq + iVqhe^"^^* C*_q 

+ i 


i SqqiC_qCq,C_qf -f 2i Tqq'CgCqiCq, . 

q' q' 


( 1 ) 


Here the first third-order term describes the interaction of 
magnon pairs {cq,c-q) and {cq',c-q'). The second third- 
order term results in a frequency shift caused by the decreas¬ 
ing of 47rMs due to the increase of the number of parametric 
magnons. The coupling between the magnon pair {cq,c-q) 
and the pumping field h is represented by the parameter Vq'. 

Vq = ^e^^‘^'>sin^eq, ( 2 ) 

4a;g 

where ujm = 747rMs, 7 = 2.8MHz/Oe, 9q = Z(q, H) and 
ipq = Z(q, q_L) are the polar and the azimuthal angles of the 
spin wave Cq{9q, ipq), respectively. According to Eqs ([T]|2]i the 
magnons with 9q = 9i = Tr(2 possess the minimal threshold 
of the parametric generation h = hth, and thus are generated 
first, exactly as it was observed in our experiment. 

Above the threshold (h > /ith), the stationary amplitude of 
the initially generated waves = c“(^, 0) = c°^q{^,TT) is 
determined by the action of the pumping field h and by the 
self-action of these waves: 




Ai 

2 


2^11 


(3) 


where is the coefficient Sqq/ for the waves of the initially 
excited parametric group Ci. Thus, the total pumping field 
in the sample consists of the external pumping field h and an 
internal pumping field induced by the magnon group Ci. This 
pumping acts not only on the waves ci but on all other waves 
in the system. If its value exceeds the excitation threshold of 
any wave with 0 ^ 7r/2 this second group C2 will be excited. 

After the excitation of the second magnon group the total 
pumping will consist of three components: The action of the 
external pumping field and the reactions of the first and the 
second groups. The effects of this complex pumping scenario 
are manifold. For example, the pumping can lead to the ex¬ 
citation of a third magnon group. It can also lead to the sup¬ 
pression of the first group of parametric magnons. The total 
pumping that drives the first group equals 

|Pl|2 = (hVlf + {SM + S2lN2f 

/ - (4) 

- 2{SiiNi + S2iN2)^J{hVi)^ - r?, 


where the product of A2 = 2|c“p and S'21 describes the ac¬ 
tion of the second group on the first one. The minus sign 
before the last term in Eq.|4]denotes the decrease of the pump 
acting on the first group. If |Pi p < [Ei p the first group will 
decay and only the second group will remain in course of time. 
This process can be understood as a stimulated scattering of 
the initially injected magnons to a new spectral area. 

In order to verify the proposed model we performed an ad¬ 
ditional experiment on detection of the second magnon group 
using the angle-resolving BLS setup shown in Fig.llla). The 
wavevector selection was realized by changing the angles 0^ 1 
and 0q|| between the sample and the incident light beam 0241] . 
It allowed us to increase the detectable wavenumber range 
up to 1.5 • 10^ rad/cm for both transversal and longitudinal 
magnons. The wavevector resolution was set to 2 • 10^ rad/cm 
by a pinhole in the path to the interferometer. To increase the 
intensity of the magnon-dependent BLS signal the separation 
between pumping pulses was decreased to 5 p,s, which was 
still sufficiently long to ensure the relaxation of the magnon 
system to the ground state after each pumping event. Accord¬ 
ingly, the duration of the pump pulses was decreased to 100 ns 
and the measurements were performed at the maximal pump 
power. 

Figure 13 a) shows the color-coded BLS intensity map mea¬ 
sured in time and wavevector space in the frequency range 
from 6.5 GHz to 7.5 GHz for magnons with q _L H. The 
two dashed lines indicate the 100 ns long microwave pump¬ 
ing pulse. The parametrically injected magnons with a fre¬ 
quency of /p /2 appear first at wavenumbers in the range from 
1.5 • 10^ rad/cm to 3 • 10"* rad/cm (see the red spot in the inten¬ 
sity map). Such a broad excitation range is due to the flatness 
of the perpendicular magnon branch near the FMR frequency 
(see Fig. [Hb), where the shadowed red segment covers the 
wavenumber range of the initially excited magnons). In the 
course of the next approximately 5 ns, some of the injected 
magnons redistribute towards smaller and higher wavenum¬ 
bers. The short scattering time is a consequence of the high 
rate of magnon-scattering processes due to the large popula¬ 
tion of the primary group of magnons. Similar to the previous 
experiment [see Fig. |2c) and (d)], the BLS signal from the 
parametrically injected magnons rapidly decreases and, ap¬ 
proximately 20 ns after the start of the pumping pulse, the 
first magnon group has completely vanished. At the same 
time, a weak BLS signal, which appears at wavenumbers be¬ 
tween 1.10 • 10^ rad/cm and 1.35 ■ 10® rad/cm and which is 
visible up to the end of the pump pulse, clearly provide ev¬ 
idence for the excitation of the second magnon group. It 
is noteworthy that the dynamics of the magnon groups ci 
and C2 [see Fig. |3d)], which was numerically calculated us¬ 
ing Eq. ([U neglecting their interaction with all other modes 
of the magnon gas, qualitatively corresponds to the experi¬ 
mental time-dependent profiles of magnon densities shown in 
Fig.|3c) as well in Figs|2jc)-(d). 

Let us address now the formation of the magnon BEC and 
specifically a question regarding the influence of the second 
magnon group on the transition of thermalized magnons to 
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FIG. 3. (color online). Color-coded logarithmic BLS intensity maps 
of (a) parametrically injected magnons (q ± H) at /p/2 and (b) the 
magnon gas (q || H) at the bottom of the spin-wave spectrum, (c) 
Change in the density of the first (red curve) and the second (green 
curve) magnon groups extracted from the intensity map in panel (a) 
for the given wavenumbers q±. (d) Time-dependent variations in 
the densities |ci|^ and |c 2 |^ of the first and the secondary groups of 
parametric magnons numerically calculated using Eq. (dJ. The long 
decay tail of the second magnon group C 2 is due to neglecting the 
nonlinear decay caused by the conventional four-magnon scattering 
processes. The intensity oscillations are caused by aperiodic dephas¬ 
ing between this magnon group and the total pumping field. 


the low-energy states. Figure [^b) represents a wavenumber- 
resolved BLS intensity map of magnons with q || H within 
the frequency range from 4.7 GHz to 5.7 GHz, which in¬ 
cludes the bottom of the magnon spectrum. One can see 
that the first thermalized magnons appear already 30 ns af¬ 
ter the rising flank of the pumping pulse. However, during 
the pump pulse the formation of the magnon BEC is pre¬ 
vented by the extremely high effective temperature of the 
parametrically driven magnon gas Hi]. After the end of the 
pumping, the magnon temperature rapidly decreases due to 
the evaporative supercooling process 11^ and magnons form 
the Bose-Einstein condensate in the global energy minimum 
at 4 ■ 10"^rad/cm [see Eig. [Tfb)]. By comparing Eig. [3a) 
and Eig. |3b) one can see that the appearance of thermal¬ 
ized magnons at the bottom of the spin-wave spectrum cor¬ 
relates with the excitation of the second group of parametric 
magnons with a perpendicular wavevector component q±^ be¬ 
yond 1.1 ■ 10® rad/cm. 

This correlation can be understood from the momentum 


conservation law. It is obvious that no magnon-magnon scat¬ 
tering process within the initially excited transversal magnon 
group is able to turn the magnon wavevector from q _L H to 
q II H orientation. At the same time, such wavevector conver¬ 
sion can easily occur if the second spin-wave group possesses 
a sufficiently large wavevector component along H. 

It is hard to predict the polar angle of the second group 
in our case as no general solution of this problem exists for 
a magnetic film of finite thickness. However, it is expected 
that this angle significantly differs from 7r/2. Eor example, 
in an unbounded medium the secondary wave group could 
be excited under the angle 6*2 — 50° 12511 . The situation is 
even more complicated for a spatially localized pumping real¬ 
ized as in our case by the narrow microstrip resonator; Be¬ 
cause excited waves can leak from the pumping area their 
losses strongly depend on their group velocities Vg^. The first 
magnon group with 0i = tt /2 propagates along the microstrip 
resonator and, thus, does not suffer any radiation losses. Prac¬ 
tically all other magnons with 9q 7^ 7r/2 can easily leak out of 
the 50 /rm wide pumping area and their excitation thresholds 
rise up significantly jloR . However, there are magnons with 
6q 7^ -n 12 which group velocities are still directed along the 
microstrip resonator. Eour of such spectral positions (n = 12, 
Oq = 35°, q = 2.16 • lO'^rad/cm, q± = ±1.2 • 10®rad/cm) 
forming the second magnon group are marked by the green 
dots on the isofrequency plane /p/2 in Eig.lHb). In general, a 
few thickness modes n can be excited simultaneously as evi¬ 
denced by the broad BLS signal produced by the second group 
in Eig. 13 a). As a result, the stimulated scattering brakes the 
orthogonality between the first magnon group {9q = 90°) and 
magnons from the bottom of the spin-wave spectrum {9q = 0), 
and thus create the prerequisites for Bose-Einstein conden¬ 
sation. Our evaluation shows that the second magnon group 
can further scatter by the conventional four-magnon scattering 
between the neighboring modes while conserving energy and 
momentum. Such intermodal scattering processes are allowed 
down to the minimum of the first n = 0 mode where the BEC 
is formed at q = 4 • 10^ rad/cm and 9q = 0 [see Eig.lTfb)]. 

In conclusion, our experimental findings suggest distinctly 
different phases of the thermalization process of a paramet¬ 
rically driven magnon gas; Eirst, energy degenerated, para¬ 
metrically magnons are generated, which, secondly, scatter 
in a multistage magnon transfer process to the bottom of 
the magnon spectrum. Thus, the parametrically stimulated 
magnon thermalization is assumed to be the essential precon¬ 
dition for the formation of the magnon Bose-Einstein conden¬ 
sate. 

We thank T. Bracher for fruitful discussions. Einancial sup¬ 
port from the Deutsche Eorschungsgemeinschaft within the 
SEB/TR 49 and from the Ukrainian Eund for Eundamental 
Research is gratefully acknowledged. 


* Clausen@physik.uni-kl.de 























5 


[1] E. Schlomann, JJ. Green, and U. Milano, J. Appl. Phys. 31, 
386S (1960). 

[2] S.M. Rezende and F.M. de Aguiar, IEEE Proc. 78, 6 (1990). 

[3] S.O. Demokritov, A.A. Serga, V.E. Demidov, B. Hillebrands, 
M.P. Kostylev, and B.A. Kalinikos, Nature 426, 159 (2003). 

[4] T. Bracher, P. Pirro, B. Obry, B. Leven, A.A. Serga, and 
B. Hillebrands, Appl. Phys. Lett. 99, 162501 (2011). 

[5] B.A. Kalinikos and M.P. Kostylev, IEEE Transactions on Mag¬ 
netics 33, 3445 (1994). 

[6] A.V. Bagada, G.A. Melkov, A.A. Serga, and A.N. Slavin, Phys. 
Rev. Lett. 79, 2137 (1997). 

[7] T. Bracher, F. Heussner, P. Pirro, T. Fischer, M. Geilen, 
B. Heinz, B. Lagel, A.A. Serga, and B. Hillebrands, Appl. Phys. 
Lett. 105, 232409 (2014). 

[8] G.A. Melkov, Yu.V. Kobljanskyj, A.A. Serga, A.N. Slavin, and 
V.S. Tiherkevich, Phys. Rev. Lett. 86, 4918 (2001). 

[9] A.V. Chumak, V.I. Vasyuchka, A.A. Serga, M.P. Kostylev, 
V.S. Tiherkevich, and B. Hillebrands, Phys. Rev. Lett. 108, 
257207 (2012). 

[10] S.O. Demokritov, V.E. Demidov, O. Dzyapko, G.A. Melkov, 
A.A. Serga, B. Hillebrands, and A.N. Slavin, Nature 443, 430 
(2006). 

[11] F. Li, W.M. Saslow, and VL. Pokrovsky, Sci. Rep. 3, 1372 
(2013). 

[12] V.E. Demidov, O. Dzyapko, S.O. Demokritov, G.A. Melkov, 
and A.N. Slavin, Phys. Rev. Lett. 99, 037205 (2007). 

[13] V.E. Demidov, O. Dzyapko, M. Buchmeier, T. Stockhoff, 
G. Schmitz, G.A. Melkov, and S.O. Demokritov, Phys. Rev. 
Lett. 101, 257201 (2008). 

[14] V. Demidov, O. Dzyapko, S. Demokritov, G. Melkov, and 


A. Slavin, Phys. Rev. Lett. 100, 047205 (2008). 

[15] A. Chumak, G. Melkov, V. Demidov, O. Dzyapko, V. Safonov, 
and S. Demokritov, Phys. Rev. Lett. 102, 187205 (2009). 

[16] A.A. Serga, V.S. Tiherkevich, C.W. Sandweg, V.I. Vasyuchka, 
D.A. Bozhko, A.V. Chumak, T. Neumann, B. Obry, 
G.A. Melkov, A.N. Slavin, and B. Hillebrands, Nat. Commun. 
5, 4452 (2014). 

[17] J. Hick, T. Kloss, and P. Kopietz, Phys. Rev. B 86, 184417 

( 2012 ). 

[18] B.A. Kalinikos and A.N. Slavin, 1. Phys. C: Solid State Phys. 
19, 7013 (1986). 

[19] A.A. Serga, C.W. Sandweg, V.I. Vasyuchka, M.B. Jungfleisch, 

B. Hillebrands, A. Kreisel, P. Kopietz, and M.P. Kostylev, Phys. 
Rev. B 86, 134403 (2012). 

[20] G.A. Melkov, V.L. Safonov, A.Y. Taranenko, and S.V Sholom, 
J. Magn. Magn. Mater. 132, 180 (1994). 

[21] G.A. Melkov and S.V Sholom, Sov. Phys. JETP 72, 341 (1991). 

[22] O. Biittner, M. Bauer, S.O. Demokritov, B. Hillebrands, 
Yu.S. Kivshar, V. Grimalsky, Yu. Rapoport, and A.N. Slavin, 
Phys. Rev. B 61, 11576 (2000). 

[23] A.G. Gurevich and G.A. Melkov, Magnetization Oscillations 
and Waves (CRC Press, New York, 1996). 

[24] C.W. Sandweg, M.B. lungfleisch, V.I. Vasyucka, A.A. Serga, 
P. Clausen, H. Schultheiss, B. Hillebrands, A. Kreisel, and 
P. Kopietz, Rev. Sci. Instrum. 81, 073902 (2010). 

[25] VV. Zautkin, V.E. Zakharov, V.S. Lvov, S.L. Musher, and S.S. 
Starobinets, Sov. Phys. JETP 35, 926 (1972). 

[26] T. Neumann, A.A. Serga, V.L Vasyuchka, and B. Hillebrands, 
Appl. Phys. Lett. 94, 192502 (2009). 



